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m ' Abstract 



Small amplitude inhomogeneous plane waves are studied as they 
propagate on the free surface of a predeformed semi-infinite body made 



to ' of Bell constrained material. The predeformation corresponds to a 



finite static pure homogeneous strain. The surface wave propagates in 

a principal direction of strain and is attenuated in another principal 

direction, orthogonal to the free surface. For these waves, the secular 

equation giving the speed of propagation is established by the method 

q . of first integrals. This equation is not the same as the secular equation 

for incompressible half-spaces, even though the Bell constraint and 

the incompressibility constraint coincide in the isotropic infinitesimal 

limit. 
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1 Introduction 

The so-called 'Bell constraint' was established empirically by James F. Bell 
in the 1980s, by conducting a great body of experiments on certain metals [TJ. 
Having subjected variety of annealed metals to different finite deformations, 
he concluded that the condition tr V = 3, where V is the left Cauchy-Green 
tensor, was always satisfied for his samples. Later, Beatty and Hayes [2] 
viewed this constraint as a purely kinematical one and used the theory of 
constrained finite elasticity to model the behaviour of the materials studied 
by Bell; it turned out that their theoretical predictions agreed remarkably 
with the experimental results (a detailed account of these advances, as well as 
references to Bell's original articles, can be found in a recent review article by 
Beatty [3].) Hence a new branch of constrained finite elasticity was born and 
a substantial literature followed, on topics such as: finite homogeneous [2] 
and nonhomogeneous [3] deformations, small deformations superimposed on 
large [5], universal relations [B] and motions [TJ, nonlinear fracture mechanics 
IE], stability analyses [HI QUI Ell E2], Cauchy stress field analysis [T3], etc. 
Small amplitude waves propagating in a finitely deformed 'Bell material' 
were examined in an unbounded media for the case of homogeneous plane 
waves J5] [TT] and in a long circular cylinder for the case of torsional waves 
[33]. The object of this paper is to establish the equation giving the speed 
of waves propagating on the plane free surface of a homogeneously deformed 
half-space made of Bell material (the secular equation). 

Elastic surface waves on solid half-spaces were first studied by Rayleigh 
[15] in the context of linear isotropic elasticity. Subsequently, the theory 
of elastic surface waves evolved in a parallel manner in the fields of lin- 
ear anisotropic elasticity [THl H7J and of finite elasticity under initial stress 
[T8"l [T9"] . Interestingly, although Rayleigh treated the case of incompress- 
ibility in his paper, surface waves propagating in internally constrained, tri- 
axially prestressed, elastic materials received only tardy attention (see for 
instance Dowhaik and Ogden [20], Chadwick [21], or Rogerson [22].) Re- 
cently this paper's author proposed a method, inspired by Mozhaev [23], to 
derive the explicit secular equation for surface acoustic waves in monoclinic 
elastic crystals [24|. In that paper, it was claimed that taking the first inte- 
grals of the tractions rather than of the displacements (as in [23]) to obtain 
the secular equation was a procedure that could 'easily accommodate internal 
constraints.' Here this claim is validated by the consideration of the Bell con- 
straint for small amplitude surface waves travelling in a principal direction 
of a homogeneously deformed half-space. 

In section 2, the incremental equations corresponding to a small ampli- 
tude displacement superposed on a large deformation [5] are recalled. These 



equations are specialized to the consideration of a surface wave in a deformed 
Bell material in section 3, where the equations of motion and the boundary 
conditions are written for the tractions on planes parallel to the free surface. 
Next in section 4, the secular equation is obtained from the vanishing of a 
2x2 certain determinant. This secular equation is compared to that ob- 
tained for incompressible finitely prestressed bodies J2U] and it is found that 
one secular equation may not be deduced from the other. However, because 
a Bell material behaves in infinitesimal motions like an incompressible one 
[2], the secular equation written when the half-space is undeformed should 
coincide with that found by Rayleigh for incompressible isotropic materials 
- it is proved that such is indeed the case. 

2 Small motion superposed on a large defor- 
mation 

Beatty and Hayes established the equations governing the elastic behaviour 
of a Bell constrained material subjected to finite deformations [2j Hj, as well 
as those describing small deformations superimposed on large |14j . Here their 
results are recalled in the context of small amplitude plane waves propagating 
in a homogeneously deformed Bell constrained half-space. 

For a general hyperelastic Bell material, the strain energy density E de- 
pends only on I 2 and I3, the respective second and third invariants of the 
left Cauchy-Green tensor V, 

£ = E(/ 2 , I 3 ), where I 2 = [(tr V) 2 - tr (V 2 )]/2, I 3 = det V. (2.1) 

Also, the Bell constraint imposes that I\, the first invariant of V, satisfy 

I t = tr V = 3, (2.2) 

for all deformations and at all times. Because of this constraint, a workless 
stress pV is introduced in the stress-strain relationship, where the scalar p 
is to be determined from the equations of motion and boundary and initial 
conditions. Hence the constitutive equation relates the Cauchy stress tensor 
T to the left Cauchy-Green tensor as [2j 

T = pY + u Q l + u 2 Y 2 . (2.3) 

Here, the material response functions uo and u 2 are defined by 

<9£ 1 as , , 



and these quantities should satisfy the Beatty-Hayes ^-inequalities [2J 

w (/2,/ 3 )<0, w 2 (/ 2 ,/ 3 )>0. (2.5) 

A solid half-space made of hyperelastic Bell material is now considered. 
In an isotropic undeformed reference state, the positions of its material par- 
ticles are denoted by X. Next, the half-space is subjected to a finite static 
pure homogeneous deformation, so that particles originally at X have moved 
to x = x(X, t). Using the principal directions of strain as a rectangular 
Cartesian coordinate system, this deformation is expressed as x a = X a X a 
(a = 1,2,3, no sum), where the principal stretch ratios Ai,A 2 ,A 3 satisfy 
Ai + A 2 + A 3 = 3, according to (j2.2p . Also, for this deformation, 

J 2 = AiA 2 + A 2 A 3 + A 3 Ai, I 3 = AiA 2 A 3 . (2-6) 

The surface x 2 = is assumed to be free of tractions and the following 
Cauchy stress tensor T satisfies the equations of equilibrium [2J, 



{T )aa = PoK + w + X 2 a u 2 , a = 1,2, 3, 

(T ) aP = 0, a ^ (3, 



(2.7) 



where u and u 2 are evaluated at I 2 , I 3 given by (12.61) . and 

Po = -(w + A^ 2 )/A 2 . (2.8) 

Next, the particle at x undergoes an incremental displacement eu(x, t) 
(with associated incremental constraint parameter ep*(x,t)), so that this 
further motion is of the form x + eu(x, t) (and the constraint parameter is of 
the form p + ep*(x, £)), where e is a small parameter. Throughout the rest 
of the paper, terms of order higher than e are neglected. The corresponding 
incremental Cauchy stress T* is given by [5]: 

T u = Aip + (AiPo + C n )- h Ci 2 - h C 13 -— , 

ax\ 0x2 0x3 

^* \ * , n dui 1 f\ l. r> \ du2 _._ n du3 
ox 1 ax 2 ax 3 

T 3 * 3 = \ 3 p* + C 31 ^ + C 32 ^ + (X 3Po + C 33 ) d ^, (2.9) 

ax\ ox 2 ox 3 

T^ &3 (A^ + A^), 

T * _i, (X 2du 2 2 du 3 
^ 2 3-&i(A 3 ^ + A 2 ^-), 

m „ , /,9<9mi ^ 9 du 3 . 



with 

Cap = 2\ a 5 a puj2 — A /3 (o;o2 + ^0^22) + AiA 2 A 3 (c<;o3 + A Q u; 2 3), ,„ ^\ 

6 7 = u 2 +Po/(3 - A 7 ), 

where the derivatives u r, w 2 r (r = 2, 3) of the material response functions 
w , u 2 are taken with respect to I? and evaluated at J 2 , I3 given by (12.61) . 
The incremental Bell constraint yields 

A^ + A^i + Ai^O; (2.11) 

OX\ OX2 OX3 

and the components of the incremental traction on planes parallel to the free 
surface X2 = are [5] 

*u = T i*2 - (AiPo + w + A^)— A 

^22 = -^22) (2.12) 

^32 = T 32 - (A3P0 + W + AJjufc)^. 

Finally, the incremental equations of motion read 

dT* d 2 u 

-dx~ =P ^ (2 ' 13) 

where p is the material density of the half-space, when it is maintained in 
the state of static homogeneous deformation. 

3 Principal surface waves 

Now the analysis of small amplitude motions in a deformed Bell half-space 
is specialized to the consideration of a plane wave propagating sinusoidally 
on the free surface X2 = in the xi -direct ion, with attenuation in the x 2 - 
direction. Without loss of generality, the component M3 is taken to be zero. 
Hence the wave may be modeled as 

u r (a;i, x 2 , t) = U T (x 2 )e lk ^- vt \ p*(x 1: x 2 , t) = kP(x 2 )e ik ^~ vt \ (3.1) 

where T = 1,2, A; is the real wave number, v is the real wave speed of 
propagation, and the specific dependence of the amplitudes Ut and P on 
x 2 need not be specified [23] • Consequently, the incremental Bell constraint 
(12. lip now reduces to 

iAitfi + A2C/2 = ° 5 ( 3 - 2 ) 



where, here and henceforward, the prime denotes differentiation with respect 
to kx 2 . The other incremental expressions recalled in the previous section 
are also simplified. 

The incremental Cauchy stress components are now given by 

T* x = k[X x P + t(XiPo + C lx )Ui + C 12 U' 2 \i 

t; 2 = k[\ 2 p + ic 21 Ut + (x 2Po + c 22 )u' 2 \, 

T 3 * 3 = k[X 3 P + iC 31 Ux + C 32 )U 2 ], (3.3) 

T* 2 = kb 3 (X 2 2 U[+iXlU 2 ), 

t; 3 = t* 3 = o. 

The third incremental equation of motion (12.131) ^ is automatically satisfied, 
while the two others reduce to 

iT* T + (T* 2 )' = -kpv 2 U r , (r = 1, 2). (3.4) 

Also, the expressions for the tractions t* 2 , t 22 , t 23 (12.121) have reduced to 

f i2 = T *2 -ik(X 1 p + uo + X\io 2 )U 2 , t* 22 = T 22 , t* 23 =T 23 = 0. (3.5) 

Now a system of four first order differential equations can be written, for 
the nonzero displacements Ui, U 2 , and for the quantities ti, t 2 , defined in 
terms of the tractions as 

t x = k-H* 12 , t 2 = krH* 22 . (3.6) 

Taking equations (12.81) . (12.101) . and (I3.2p into account, this system can be 
written as 

t[ + iXiX 2 l t 2 - (X^C - pv 2 )XJ x = 0, 

t 2 + ih - [6 3 (A? - Xl) - pv 2 }U 2 = 0, (3.7) 

U^ + iX 1 X 2 1 U 1 = 0, 

b 3 X 2 2 U[ + ib 3 X 2 2 U 2 - ti = 0, 

where C is defined by 

C = A 1 A2C11 + AiA 2 C22 ~~ W2 ~~ ^21 ~~ 2cjg — [X 1 + X 2 )cj 2 , (3-8) 

and the following expression for 6 3 A|, obtained from (12.101) . and (12.81) . has 
also been used, 

b 3 X 2 2 = -{u - X l X 2 u 2 )/(l + A^ 1 ). (3.9) 



The system must be solved when subject to the following boundary con- 
ditions. First, the surface x 2 = remains free of traction, so that 

t* 2 (0)=0, (r = l,2); (3.10) 

and second, the displacement must vanish at infinite distance from the free 
surface, so that 

U r (oo) = 0, P(oo)=0, (r = l,2). (3.11) 

4 Secular equation 

Here the secular equation is derived by applying the method of first integrals 
to the system (13. 7p . 

Proceed as follows. First, rewrite the boundary conditions ( 13 . 1 [) - ( 13 . 1 1 [) 

as 

£ r (0) = £ r (oo) = 0, (r = l,2). (4.1) 

These last conditions, together with (j3.7j) i ■?. imply that 

t' r (oo)=0, (r = l,2). (4.2) 

Next, obtain Up in terms of tp, t' r from ( I3.7p -i o. and U r in terms of Up, tp 
from (I3.7p ^/i. or equivalently, in terms of tr, tp- Finally, substitute these 
expressions into the derivatives of equations fl3.7[) i _q to get, 

[6 3 (A? - \l) - P v 2 ]t'[ + i/3 12 t' 2 - -L^A^C - pv 2 )(b 3 \l - pv 2 )t x = 0, 

3 A 2 

(4.3) 

(AxA^C - pv 2 )t'i + z/3 12 t; - \l\ 2 2 [h{\l - A!) - pv% = 0, 

where (13. 9p has been used, and /3i 2 is defined by 

/012 = A 1 A 2 - 1 [6 3 (A? - \\) + C] - (1 + A1AJV 2 . (4.4) 

The quantity (5\2 does not play any role in the secular equation, as is now 
seen. Multiply (I4.3p 1 by t[ and (14.3j) -> by t' 2 , and integrate with respect to 
kx 2 between x 2 = and x 2 = 00, to obtain, using ( 14. ip and (14. 2p . 

[b 3 (Xl - \ 2 2 ) - pv 2 }t[(0) 2 - 2zf3 l2 jt' l t> 2 = 0, 

(AiA^ 1 ^ - P v 2 )t' 2 (0) 2 - 2if3 12 ft[t' 2 = 0, (4.5) 



and therefore, 

[bs(\l - Af) - pv 2 ]t'M 2 - (AiA^C - pv 2 )t'M 2 = 0. (4.6) 

Similarly, multiply fO]b by [6 3 (A? - A^) - pw 2 ]^ + «/3i 2 £ 2 and (|4~3]) 2 by 
(AiA^C — pv 2 )t' 2 + «/3i2^i, and integrate with respect to kx 2 between x 2 = 
and x 2 = oo, to obtain, using (14. ip and (14.21) . 

[6 3 (A 2 - A 2 ) - P v 2 ]H[(0) 2 - 2ip l *(\ 1 \; 1 C- F v 2 )(b 3 \ 2 1 -pv 2 )ft 1 t 2 = 0, 

3 A 2 

(X 1 X 2 1 C - P v 2 )%(0) 2 - 2^ 12 A 2 A 2 - 2 [6 3 (A 2 - A 2 ) - pv 2 \jt x t 2 = 0, (4.7) 
and therefore, 

[^-A^-p^A^O) 2 

" 7TT2 ( WC - P^ 2 ) 3 (M 2 - P^ 2 )4(0) 2 = 0. (4.8) 
O3A2 

The two homogeneous equations (14.61) and (14. 8 p for the two unknowns 
£i(0) 2 and t' 2 (0) 2 yield nontrivial solutions only when the corresponding de- 
terminant is zero. This condition provides the secular equation for surface 
waves propagating in a principal direction of a deformed Bell material as, 

(6 3 A 2 - pv^X^C - pv 2 ) 2 = hX\[h{X\ - A 2 ) - pv 2 } 2 . (4.9) 

This equation is distinct from that obtained for deformed incompressible 
materials by Dowaikh and Ogden [20], which reads as 

(6 3 A 2 - pv 2 ){C - pv 2 ) 2 = M 2 [6 3 (A 2 - A 2 ) - pv 2 } 2 , (4.10) 

where 63 and C are defined in terms of the stretch ratios Ai, A 2 and the stored 
energy density W of the incompressible material as 

~ ow dw 2 2 

& 3-(Ai^- A 2 — )/(X 1 A 2 ), 

~ ^d 2 W , , d 2 W ^d 2 W BW . 

G = A; _- 2AlA2 __ + Ai _ + 2A2 _ (4.11, 

Although equations (14. 9 p and (I4.10p share a similar left hand-side, their right 
hand-sides are different and could not have been a priori deduced one from 
another. In particular, 63 and AiA 2 ~ 1 C can formally be transformed into 63 
and C respectively to get from the left hand-side of (14. 9p to the left hand- 
side of (14.101) . but this transformation does not work for the right hand-sides 

8 



of the equations, because the factor A^ is replaced by A|. Hence, when the 
influence of the prestrain on the propagation of surface waves is studied 
[20] , the behaviour of Bell materials is in contrast to that of incompressible 
materials, except in the special degenerate case of a bi-axial prestrain such 
that Ai = A2. 

However, when the Bell material is unstrained for the problem at hand 
(Ai = A 2 = A 3 = 1), only infinitesimal deformations (the surface waves) 
occur in the now isotropic half-space. In such a case, Beatty and Hayes [2] 
show that the Bell constraint is equivalent to the incompressibility constraint. 
Here, the quantities C, 63 are then given by 

AiA^C = C = 2(oo 2 - u ) = 46 3 , (4.12) 

and the secular equation (14. 9 p reduces to 

(l-0(4-0 2 = e 2 , or e-$e + 24^-16 = 0, (4.13) 

where £ = 2pv 2 /(u 2 —ou ). This last equation is that established by Rayleigh 
for surface waves in isotropic linearly elastic incompressible half-spaces. 



References 

[1] J.F. Bell, Contemporary perspectives in finite strain plasticity. Int. J. 
Plasticity 1, 3-27 (1985). 

[2] M.F. Beatty and M.A. Hayes, Deformations of an elastic, internally 
constrained material. Part 1: Homogeneous deformations. J. Elasticity 
29, 1-84 (1992). 

[3] M.F. Beatty, Hyperelastic Bell materials: retrospection, experiments, 
theory, In: Nonlinear Elasticity: Theory and Applications. Eds: Y.B. Fu 
and R.W. Ogden, 58-134 (Cambridge University Press, London 2001). 

[4] M.F. Beatty and M.A. Hayes, Deformations of an elastic, internally 
constrained material. Part 2: Non-homogeneous deformations. Q. Jl. 
Mech. Appl. Math. 45, 663-709 (1992). 

[5] M.F. Beatty and M.A. Hayes, Deformations of an elastic, internally con- 
strained material. Part 3: Small superimposed deformations and waves. 
Z. angew. Math. Phys. 46, 72-106 (1995). 

[6] E. Pucci and G. Saccomandi, Universal relations in constrained elastic- 
ity. Math. Mech. Solids 1, 207-217 (1996). 

9 



[7] E. Pucci and G. Saccomandi, Universal motions for constrained simple 
materials. Int. J. N on- Linear Mech. 34, 469-484 (1999). 

[8] A.M. Tarantino, Nonlinear fracture mechanics for an elastic Bell mate- 
rial. Q. Jl. Mech. Appl. Math. 50, 435-456 (1997). 

[9] F.X. Pan and M.F. Beatty, Instability of a Bell constrained cylindrical 
tube under end thrust. Part 1: Theoretical development. Math. Mech. 
Solids 2, 243-273 (1997). 

[10] M.F. Beatty and F.X. Pan, Stability of an internally constrained, hy- 
perelastic slab. Int. J. Non-Linear Mech. 33, 867-906 (1998). 

[11] F.X. Pan and M.F. Beatty, Instability of an internally constrained hy- 
perelastic material. Int. J. Non-Linear Mech. 34, 169-177 (1999). 

[12] F.X. Pan and M.F. Beatty, Instability of a Bell constrained cylindrical 
tube under end thrust. Part 2: Examples, thin tube analysis. Math. 
Mech. Solids 4, 227-250 (1999). 

[13] M. Hayes and G. Saccomandi, The Cauchy stress tensor for a material 
subject to an isotropic internal constraint. J. Eng. Math. 37, 85-92 
(2000). 

[14] M.F. Beatty and M.A. Hayes, Small amplitude torsional waves prop- 
agating in a Bell material, In: Nonlinear Waves in Solids. Eds: J.L. 
Wegner and F.R. Norwood, 67-72 (ASME 1995). 

[15] Lord Rayleigh, On waves propagated along the plane surface of an elastic 
solid. Proc. R. Soc. London A17, 4-11 (1885). 

[16] A.N. Stroh, Dislocations and cracks in anisotropic elasticity. Phil. Mag. 
3, 625-646 (1958). 

[17] A.N. Stroh, Steady state problems in anisotropic elasticity. J. Math. 
Phys. 41, 77-103 (1962). 

[18] M.A. Biot, The influence of initial stress on elastic waves. J. Appl. Phys. 
11, 522-530 (1940). 

[19] M. Hayes and R.S. Rivlin, Surface waves in deformed elastic materials. 
Arch. Rati. Mech. Anal. 8, 358-380 (1961). 

[20] M.A. Dowaikh and R.W. Ogden, On surface waves and deformations 
in a pre-stressed incompressible elastic solid. IMA J. Appl. Math. 44, 
261-284 (1990). 

10 



[21] P. Chadwick, The application of the Stroh formalism to prestressed elas- 
tic media. Math. Mech. Solids 97, 379-403 (1997). 

[22] G.A. Rogerson, On the existence of surface waves and the propagation of 
plate waves in pre-stressed fibre-reinforced composites. J. Mech. Phys. 
Solids 46, 1581-1612 (1998). 

[23] V.G. Mozhaev, Some new ideas in the theory of surface acoustic waves 
in anisotropic media, In: IUTAM Symposium on anisotropy, inhomo- 
geneity and nonlinearity in solids. Eds: D.F. Parker and A.H. England, 
455-462 (Kluwer, Holland 1994). 

[24] M. Destrade, The explicit secular equation for surface acoustic waves in 
monoclinic elastic crystals. J. Acoust. Soc. Am. 109, 1398-1402 (2001). 

[25] S. Nair and D.A. Sotiropoulos, Elastic waves in orthotropic incompress- 
ible materials and reflection from an interface. J. Acoust. Soc. Am. 102, 
102-109 (1997). 

[26] D.A. Sotiropoulos and S. Nair, Elastic waves in monoclinic incompress- 
ible materials and reflection from an interface, J. Acoust. Soc. Am. 105, 
2981-2983 (1999). 

[27] M.Y. Yu, Surface polaritons in nonlinear media, Phys. Rev. A28, 1855- 
1856 (1987). 



11 



